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The equations-of-motion for the density matrix are derived in a multiband model to describe the 
response of semiconductors (bulk or quantum well structures) under optical excitation with arbi- 
trary polarization. The multiband model used, comprising the twofold conduction band and the 
fourfold topmost valence band (or heavy- and light-hole states), incorporates spin-splitting of the 
single-particle states. The interaction terms include besides the direct Coulomb coupling between 
carriers also the electron-hole exchange interaction, which together with the spin-splitting terms 
is responsible for spin relaxation. Applying the Hartree-Fock truncation scheme leads to a set of 
coherent semiconductor Bloch equations for the multiband case. This concept provides the theoreti- 
cal frame for describing phenomena connected with optical response under excitation with arbitrary 
light polarization and spin relaxation: polarized optical response, polarization dynamics of VCSELs, 
spin relaxation, and the circular photovoltaic effect. 



PACS Num. 72.25.Fe, 72.25.Rb, 78.47.+p 
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I. INTRODUCTION 



The spin degree-of- freedom of electrons has recently attracted increasing attention in the context of spintronics. 1 
In semiconductor quantum structures it is intimately connected with circularly polarized light that is frequently 
used to create spin-polarized carriers and to detect their lifetime with respect to spin relaxation in pump-probe 
experiments. 2-4 The investigation of optical orientation and of the mechanisms of spin relaxation, first carried out 
in bulk semiconductors 5 , have become the basis for ongoing studies in semiconductor quantum wells (QW), which 
are the likely material structure to realize the spin transistor. 6 Theoretical work devoted to spin relaxation makes 
use preferentially of the density matrix to describe the thermal expectation value of the spin as an observable. 7 For 
spin-polarized electrons in the conduction band of a bulk semiconductor or in the lowest electron subband of a QW 
structure it is sufficient to have the 2x2 spin-density matrix and its time evolution due to generation/annihilation and 
scattering processes. However, an adequate modelling of polarized photoluminesccncc involving intcrband transitions, 
i.e. the calculation of the thermal expectation value of the dipole operator for circularly polarized light, requires the 
knowledge of the density matrix for conduction- and valence-band states in a multiband scheme. 8 The dependence on 
the intensity of the exciting light is an important experimental aspect and should be included in the density matrix. 9 

For vertical-cavity surface-emitting lasers (VCSEL) the phenomenon of polarization instability has recently become 
a topic of interest. 10,11 The VCSEL geometry allows for two independent polarizations of the laser light, while VCSEL 
modelling, using the optical or semiconductor Bloch equations (SBE), was based so far on the two- level system that 
does not account for polarization degrees-of-freedom. 12,13 This deficiency has been removed partially by formulating a 
phenomcnological set of optical (or Maxwell) Bloch equations for a pair of two-level systems, that can emit separately 
light with right or left circular polarization and which are coupled by spin-flip processes. 10 

Semiconductor Bloch equations can be studied in different time regimes: the coherent regime - when the interband 
polarization is in phase with the exciting electromagnetic field - is characterized by the Rabi oscillations, while the 
quasi-equilibrium situation is described by the optical (linear or nonlinear) susceptibility. In extending the two-level 
system, showing one Rabi frequency, to a three (or actually six) level system, comprising the spin-degenerate conduc- 
tion, heavy- and light-hole bands, Binder and Lindbcrg 14 formulated a set of equations for the optical polarization 
functions and carrier distributions to describe coherent nonlinear response of this more realistic quantum- well model 
under excitation with circularly polarized light. This treatment includes band mixing due to heavy- light hole cou- 
pling in the Luttinger Hamiltonian and many-body effects in a Hartree-Fock(HF) truncation with single-particle self 
energies determined by the direct Coulomb interaction. This concept has been applied successfully in studies of the 
intervalence band coherence due to coupled optical Stark shifts. 15 However, spin-flip processes that become possible 
due to band structure effects or due to electron-hole exchange (not in the scope of these studies) were not considered. 
Several more recent papers reported on the observation of the circular photogalvanic effect (CPGE) in n- and p-doped 
semiconductor quantum well structures. 16 This effect, originally predicted for bulk semiconductors 17 , is based on con- 
verting the helicity of light into a directed stationary current. It represents a nonlinear response of the system to the 
intense exciting light which can be described by a third rank tensor. The effect is caused by creating a nonequilibrium 
carrier distribution in the spin-split subbands by excitation with circularly polarized light. The dependence of the 
effect on the light intensity, i.e. its saturation behavior, bears information on the spin-relaxation and should allow 
to extract the characteristic times. 18 A microscopic description of these phenomena requires the formulation of the 
nonlinear response under excitation with circular (in general elliptic) polarization. 

Among the different mechanisms of spin relaxation 19 those named after D'yakonov and Perel' (DP) 7 and Bir, 
Aronov, and Pikus (BAP) 20 seem to dominate in semiconductor structures not containing magnetic impurities. The 
DP mechanism is a single-particle effect. It is intimately connected with spin-splitting of the electronic states caused 
by spin-orbit coupling and inversion asymmetry. The BAP mechanism results from electron-hole exchange scattering 
and depends on the carrier concentration resulting from doping or intense optical excitation. Both mechanisms are of 
the motional-narrowing type, i.e. the spin-relaxation rate is proportional to the momentum scattering time. 21,22 In 
experiments with circularly polarized light in the visible spectral range electron-hole pairs are created (bipolar optical 
orientation) and spin relaxation takes place due to both mechanisms, a situation which was not accounted for so far 
in theoretical work. The BAP mechanism can be avoided in far-infrared experiments with monopolar spin orientation 
connected with transitions between electron (or hole) subbands. 16,18 

In order to provide a general theoretical frame for all these phenomena the equations-of-motion for the density 
matrix are formulated in a multiband model with the spin and polarization degrees-of-frccdom being taken into 
account. The system Hamiltonian 

H = H + H'(t)+H Co ui (1) 

consists of the single-particle part, H + H'(t), describing the involved electron states and the interaction with the 
electromagnetic field, and the Coulomb interaction Hcoui- The multiband model comprises the twofold conduction, 
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heavy- and light-hole bands (or the respective lowest subbands in a QW structure) as in Ref. 14 but includes the 
mechanisms of spin-splitting as origin of the DP mechanism. In Section II the equations-of-motion will be formulated 
by considering only the single-particle part of H in order to introduce the multiband model and to get familiar with 
the notation. Many-body effects (due to free carriers introduced by doping or created by optical excitation) are 
considered in Section III. Besides the direct Coulomb interaction between carriers also the electron-hole exchange 
coupling is taken into account. It gives rise to spin- flip processes and represents the origin of the BAP mechanism 
of spin- relaxation. By following Refs. 12-14 the Hartree-Fock truncation is applied in Section IV to get a closed set 
of equations for the elements of the density matrix that describe the coherent part of the problem. They turn out 
to be more general than those known from the literature 8 ' 14 by taking into account both the polarization and spin 
degrees-of-freedom. In Section V these equations will be discussed with respect to simplifications (and including 
phenomenological damping) which allow to recover published results and to point out applications to new nonlinear 
response phenomena addressed here. 

II. THE SINGLE-PARTICLE PART 

The single-particle Hamiltonian H is conveniently written in the second quantized form 

Ho=J2 ^Jk)ai ma (k)a ama (k) (2) 

am a k 

where a,m a ,k are the quantum numbers of a complete set of eigenstates with energy e°J^ and a 1 * (a) are cre- 
ation(annihilation) operators of the corresponding states. For bulk semiconductors a denotes the energy bands, m a 
the degeneracy within a given band, and k the three-dimensional (3D) wave vector. In the following we restrict the 
basis to the lowest conduction band (a = c) with m c = ±1/2 and the topmost valence band (a — v) with m v = ±3/2 
or ±1/2 for the heavy and light holes, respectively. The notation applies as well to QW structures with the lowest 
electron (a — c,m c = ±1/2), heavy (a = v, m v = ±3/2), and light hole (a = v, m v = ±1/2) subbands and a 2D wave 
vector k. (An extension to more subbands would require book keeping of subband indices, which is avoided here.) The 
single-particle energies (k) are obtained by diagonalizing a 2 x 2 (4 x 4) Hamiltonian for the conduction (valence) 
band states, which may include besides the diagonal free-particle kinetic energy also nonparabolic corrections from 
a higher order k -p decoupling. 23 This includes terms resulting from bulk-inversion asymmetry or from the asymme- 
try of the quantum well and lead to a removal of the spin-degeneracy of the (sub)bands. 24 Thus, the single-particle 
wave functions are expanded in a basis of band-edge (k = 0) Bloch functions u am > a (x) (x comprising space and spin 
coordinates) which for the bulk semiconductor reads 

* ama k(30 = e rk -^C« am ,Jk)u arn , a (x) . (3) 

The coefficients C^ Qm / (k) describe within a band a the mixing of the band-edge states at finite k; for the valence 
band this includes coupling between heavy- and light-hole states. Note, that the angular-momentum classification 
is exact only at k = 0, where C^ am , (0) = 5 mam > a . For QW structures (with z being the growth direction) the 
single-particle wave function takes the form 

V> aroaj g(z) = e^^C ( ^;\k,z)u am , a (x) . (4) 

Here k — (k x ,k y ,0) is the in-plane wave vector, f = (x,y,0), and C^ a ^{k, z) are the subband functions including k 

dependent mixing, while for k — > we have C^°^{k, z) — > Cin/°^(0; z )^m a m' a - 

With this formulation of the single-particle states Hq is block-diagonal in the band indices a = c, v and it is 
convenient to replace a C m a {k) — > c mc (k) for electrons and a vmv (k) — > v^_ m (— k) for holes (adopting the picture of 
holes as time-reversed electron states) and similar for the hermitian adjoint fermion operators. In this notation (using 
the Kramers degeneracy e" ma (— k) = £? 1la {k)) H (Eq. (2)) reads 

#o - E £ ™. W c ™. W c »- $ + E $)v m „ (k)vl v (k) . (5) 

ro c k m v k 

Likewise the interaction with the electromagnetic field can be written 

#'(*) = - E {E(t)-dZ m Mcl c (k)vl mv (-k) + hx.} (6) 

m c m v k 
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with the matrix of the dipole operator 



in 'rn/„ 



for the bulk case (here C^, mv (k) denote the expansion coefficients of the hole wave functions) and 

dt c mM = e £ J Clnt ] *(k, z)R m , m ,(^\k, z)dz . (8) 

m'm' 



for the QW case. Here the real electric field vector is E(t) = E{uj)e lut + E*{u))e~ iujt (E*(oj) = E(-co)) and R m ' c m< v , 
the matrix element of r, is taken between the band-edge Bloch functions u cm ' c (x) and u vm ' v (x) (which are angular- 
momentum eigenfunctions) . For light with arbitrary polarization propagating in the z-direction the complex field 
amplitudes can be written E±(ui) = E (u))(e x + e y e ±4 ^)/2 where the ±-sign indicates (for <fi ^ 0) the hclicity of light. 
For 4> = it/2 this is the right and left circular polarization. The matrix elements of E±(uj) ■ R m 'm' are given in the 
following table (R =< S\x\X > /y/2): 

m c = +h -^(l + ie^*) fa£(l -ie ±irt> ) 
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m v = 



2 V- ' ) " 2V3 



_|n*(l + ie ±<*) ■ ^«(l_ ie ±i*) 



i — i n nr 

2 2 2 2 



The single-particle Hamiltonian H + H'(t) is a sum of contributions from different wave vectors k at each of which 
we have a six-level system according to the multiband model. In analogy with the two-level model of Refs. 12,13 one 
can now formulate the equations-of-motion for the operators Cfh c ,V- mv v^_ m , and v- mv Cm c 



(ihd t - e c mc (k) + e% lc (k)) cl e (k, t)c mc (k, t) = 

= £ {E(t) ■ dZ m Mcl c (k,t)vl mv (-k,t) - h.c) (9) 

m v 

(ihd t +e v _ mv {-k) - el mv (-k))v. mv (-k,t)vU v (-k,t) = 

{E(t) ■ dZ m Mcl c (k,t)vl mv (-k,t) h.c] (10) 



E 



; E(t) ■ d^ m , (k)v^ (-k, t)vl m , (-k, i) - £ E(t) ■ d% cmv (k)cl, (k, t)Crn c (k,t). (11) 



E- 



These operators can be identified according to 

Pm c m c (M) = c L c (k,t)cfh c {k,t) , m c , rh c = ±1/2 (12) 

P™„m„(M) = v- mv (-k,i)vl mv (-k,t) , m v ,m v = ±1/2, ±3/2 (13) 

PLfn c (k,t) = V- mv (-k,t)Crn c {k,t) (14) 

with the operators p c (1 — p v ) of the spin density matrix for electrons (holes) and P of the interband polarization 
matrix. The matrix E(t) ■ <i™ cTOij (k) is the multiband generalization to the (unrenormalized) Rabi frequency of the 
two-level model. When taking the thermal average of the matrix operators Eqs. (9)-(ll) correspond to the single- 
particle parts of Eqs. (l)-(3) of Ref. 14 , which however do not include spin splitting of the electron and hole eigenstates. 
(Note, that in the present case the conduction- and valence-band Hamiltonians are diagonalized, see Eq. (5)). 

III. THE MANY-BODY TERMS 

In correspondence with the formulation of the single-particle part of the system Hamiltonian the Coulomb interaction 

Hcoui = lf d * J dx'*^x)&(x!)v(\f-7<\)*{x')*{x) (15) 
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is decomposed by splitting the field operator 

tf (z) = * e (ss) + (16) 

into its electron and hole parts, where x stands again for space and spin variables. Thus one obtains 

Hcoui = H ee + H hh + H^ h + H* h (17) 

where the individual terms have the same form for the bulk and QW case (time-dependence of the creation and 
annihilation operators is understood) 



H ee = 



\ E E V ™c 1 mc 2 mc3me 4 (M^>L 1 (£+^ ( 18 ) 

H hh = \ E E V ^ im „ 2 ^ 3m „ 4 & k', q)v mvi (k - q)v mv2 (k' + fiv^ (fc')^„ 4 (k) , (19) 

= - E E Cf^™^ & ?14 C1 (fc + ?1c mc2 (fc)^ (fc')^ (# + 9) , (20) 

^ = E E V ^m vl m C2 m V2 (k, ?")4 01 + 9>mc 2 + ?>^ 2 (fe)« m , 1 (#) . (21) 



™ci m c 2 fefe' 



They describe the interaction between electrons in the conduction band (H ee ), between holes in the valence band 
(Hhh) 7 and the direct electron-hole Coulomb interaction {H^ h ). These terms are known from previous treatments 
of the problem. 9,12 ' 14 The electron-hole exchange interaction (H* h ), however, has not been considered so far in the 
context of SBE. 25 Together with the spin-splitting mechanism in the single-particle part (Section II) it is the essential 
extension of the existing theory and will be discussed in the course of this paper. 

Due to the diagonalization of the (sub)band Hamiltonians for electrons and holes (see Eqs. (2) and (4)) the inter- 
action matrixelements depend on the expansion in terms of the band-edge Bloch functions. For the bulk system they 
read 

V e m \ ...m Ci (k, k>, q) = v(q) ]T CZ cimc $ + ?1^* C2 < (*' - O)^^, (k')C c mc4mc (k) , (22) 

m c m' c 

V% i ... m Jk,k',q)=v(q) Y ^ im „(£-^- 2TO ,^ (23) 

m v m' v 

Vt°.. mc2 {k,k',q)=v{q) Y, C^m^ + flC^^ (24) 

m c m v 

V e m h ;X. mv2 (k, k>, q) = v{q) Y, M 2 C Xl C ™ cl m c (~k + q)C^ m , (-P + q) x 

m c m v m' c m' v 

X ^m vl m v {k / )Cm V2 m' v {k) (25) 

with the expansion coefficients C c for electron and C v for hole wavefunctions. Here v(q) is the Fourier transformed 
3D Coulomb potential. For the QW case one has instead 



VZ imc2 m C3 mJk,k',q) = v(q) Y J dz J dz'e^^ 

m c m' 



X 



*C^\k + q,z)C;* ) \k>-q,z>)C;> \k',z')^\k,z) (26) 

and corresponding expressions for V /lh , V eh,c , and \? eh > x t with the band-mixing at finite in-plane k considered in the 
subband functions. Here v(q) is the 2D Coulomb potential (see e.g. 9 Eq. (A. 18)). The Matrix M in the exchange 
interaction is due to the coupling of electron and hole angular momenta, known from studies of the finestructure 
splitting of excitons 26-28 or of the BAP spin-relaxation mechanism 21,22 . Without mixing of the band-edge states 
the interaction matrixelements y ee ,V /l ' 1 , anc l yeh.c ta k e mucn simpler forms as e.g. V^ c mv mv TOc (k,k',q) — ► 

v{q)5 mcimc2 5 mviTnv2 and V^^^^ -> v(q)M™l\m? 2 for the 3D (or bulk) case and similar for the QW case. 
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By calculating the commutators of the operators from Eqs. (12)-(14) with Hcoui one finds the following additional 
terms which are to be added to the Eqs. (9)-(ll): 



Hcoui , c4 c (fc)cm c (fc)J = 



kqm' c 



(k, k, q)cl lci (k + q)c^ (k - q) CThr (k) $)) 
E { V t;Z vimv2 m'fi ~ Z k, q)cl c (k) Cm ,(k q) 



yeh.C 



' 7n' c m vi m V2 m c 

{k, k, q)& m , c (k + q)c mc (k)J (k)v mvi (k + q) 

+ E ( V ^cm vl m' c m V2 (-k + Q, k^fLMCrn'S-k + ^V^i-k + fiVrn^ik) 

- <l imcm „ 2 (*>"* + <?, Q)cl,(-k + q)c mc (k)vl V2 (k)v mvi (-k + q))}, 

= - E { E ( V ^ -rh v m> v m v , $,-%- q, q)V m ^ (k - q)V. mv (~k - (k) 

- V m vl m' v -m v m V4 ~k, g>m n (fc ~ q)V m > v (~k + qjvL^ (-£)*4„ 4 <>)) 

E (<f<-„ 2 (*> ~k, Q)v m , (-k + q)vl mv (~k) 



- V, 



1 "-C2 



m cl —m v m'm 



' m,„ (k, -k - q, q)v- mv (~k)vL (-k - q)) S (k + q)c m ^ (k) 



- V, 



E ( V m^m' v m C2 -m v (~k, K q)^ (k + q)C m ^ (-k + q)V n< (k)vl mv (~k) 

, (k, -k, q)c\ 1lci (-% + q)c mc2 (k + q)v- mv (-fc)^ (fc)) } , 



1 ""C2 



m cl —m v rn r ,. in' 



Hcoui , v- m „(-£)c™ c (£) 
= - V V ee 

/ ^ y rn ci 



m c m C o m c 



Ic l " Ic 3 c 4 



(fc, fc + q, q)v- mv (-fc)4„ ci (fc + q)c mci (k)c mc3 (k + q) 

E V ™ V2 m vl m V3 -m v (-k> -3)"m ei - «>m» 2 + $)t4„ 3 (%m t (fc) 

TTItj^ tti V 2 

fcq 

+ E {T, V tZ vimim Jk-qXqlv- m A-k)c mc2 (k 

m C2 m Vl m f^ 
kg 

+ E V <m 01 -m„m C2 (*> (* + ^ C »^2 + ^ C ™c (*)) 

+ E ( E V tt vl m C2 m' v (~k + 1, k, q)v. mv (-k)c mc2 (-% + q)v] n , (-k + q)v mvi (k) 
+ E ™ C2 -m v (~kX q)cl, (k + q)c mc2 {-% + q)v mvi (k) Cfhc {%)) . 



(27) 



(28) 



(29) 



Equations (27) and (29) reduce to Eqs. (A. 29) and (A. 28), respectively, of Ref. 9 when the band mixing and the 
exchange terms are neglected. The characteristic effect of the Coulomb interaction is to add four-operator terms to 



6 



the equations of motion of the two-operator terms, thus giving rise to a hierarchy of such equations which can be 
solved only approximately. The choice of the method will depend on the scenario that is to be described: the relevant 
time scale (to distinguish between the coherent or quasi-equilibrium regime) , the carrier density, or the power of the 
driving electric field. 12 ' 13,31 . 

IV. THE COHERENT DENSITY-MATRIX EQUATIONS 

A set of equations describing the coherent regime can be derived by following Refs. 12 ' 13 and applying the HF 
truncation to factorize the four-operator terms. By taking thermal expectation values and keeping only factors, which 
are diagonal in the wave vector 

c ln c (^it) c m c (k,t) — > (cl nc (k,t)Cfh c (k,t)) = P mr m r (k,t) 
v -rh v (-£> t)V- mv {-k, t) -» (ul m „ (-k, t)v_ mv (-k, t)} = p^ mnj (k, t) 

V_ mv {-%, t)Crn c (k, t) -» {V- mv {-k, t) Crhc {k, t)) = (k, t) (30) 

one finds a closed set of equations for the components of the density matrix 
[ifidt - e c mc {%) + e% e (fc)) p e mcrhc (k, t) = 

= + E |^*n c m„ (&J t)Pm c m v (k, t) — £l mr m v (k, t)Pfh a m v (&i 0} 

+ E { S -c< (*, *)P<m c (*, *) - P e m c m' c (*, *)S^* e (*, t) } (31) 

(ihd t (k) + (k))p h ^ mv (k, t) = 

- E { S -.K *) - ^:^„ (fc, *) } (32) 



+ 4, c (fc) - c„ (*)) PLm v (k, t) = 

mi m£ 

_ E ^m' c m^ t)^m' r m c 0, t) + E ^m„mi t)Pfh c rn' v *) ■ (33) 
m' c mi 

These equations are generalizations of Eqs. (l)-(3) in Ref. 14 in that they include the spin-splitting mechanisms in 
the single-particle part and by considering also the electron-hole exchange interaction in the electron and hole self 
energies: 

^m c m^ (£> t) — ~ E { E ^m' c fh c m c fh' c — Q)Pfii c m' c — 9i *) 

q rh c m' c 

- E <-m„m c -mi(-^+9,-fc + g,9>^ ro i(fc-9,i)} (34) 

S^„ m / (M)=-E { E V -mi-m„-m„-mi M + <?, S) (*A„mi - Pr<m„ (* - 5, *)) 

q rh v rh' v 

+ E V m c -mimi-m„ ("*> ?>m c mi + £ *) } ■ (35) 
m c m f c 
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The renormalizcd dipole matrix reads 

q,-k,q)P r * Km , v (k-q,t) (36) 



and differs from the corresponding expression in Ref. 14 (Eq. (6)) only by the alternative choice of the basis. 

The self-energy terms have to be evaluated by summing over pairs of angular momenta for electrons and holes and 
thus couple to all components of the electron or hole density-matrix. Take, e.g., the equation of motion for p% riaifLa (k, t): 
the diagonal terms with E^ cmc and S^ cJRc renormalize the single-particle energies, while the off-diagonal terms 
correspond to spin-flip scattering. The latter can be due to electron-electron interaction modified by band-mixing 
but also to electron-hole exchange interaction. In fact, by neglecting in the self-energy, Eq. (34), the band-mixing 
it can be seen that the electron-electron interaction (because of V^f, m c m a m> ~ dm' c m> c fim a m c ) does not contribute to 
spin flips while the second term (from electron- hole exchange interaction) retains the sum over m, , m' v and couples 
electron and hole angular momenta. As it flips an electron spin simultaneous with a hole spin (if holes are present), 
it can be identified as the origin of the BAP mechanism. The self-energy corrections due to V ee and V hh contain the 
spin-scattering processes mediated by band-mixing in combination with spin-orbit coupling. They are more efficient 
for holes (due to their p-character), while spin-orbit coupling in the electron states (sometimes called Rashba 29 and 
Dresselhaus term 30 ) results from band-mixing at finite k. Thus Eqs. (31)-(33) represent an extension of the coherent 
SBE 12 ' 13 to the multiband case and arbitrary polarization of the driving electromagnetic field, which goes beyond 
Ref. 1 by including the spin-flip processes. These equations provide the theoretical background to describe experiments 
in the coherent regime that address the spin-degree of freedom like spin-echo or four-wave mixing experiments with 
elliptically polarized light (in the latter case one could go beyond the HF truncation by using the dynamically controlled 
truncation scheme 31 ). 

In closing this section a compact form of the Eqs. (31)-(33) is presented which becomes possible by making use of 
matrix notation. 32 In the sixfold space of the multiband model used here the Hamiltonian Ho(fc) + S(fc) (without the 
coupling to the electromagnetic field) is block diagonal with 2x2 and 4x4 diagonal blocks for the single particle 
and self-energy parts of conduction (Hg(fc) + S c (fc)) and valence band states (Hg(fc) + £"(&)), respectively. The 

interaction with the electromagnetic field, described by the matrix — 0(fc, t) of the (renormalized) dipole interaction, 
couples between the valence and conduction band states and has off-diagonal blocks only. Similarly, the 6x6 density 
matrix p(k, t) consists of diagonal blocks, the 2x2 electron density matrix p e and the 4x4 hole density matrix 1— p h , 
and off-diagonal blocks for the interband polarization P (the upper right 2x4 block) and its hermitian conjugate. 
Written with these matrices Eqs. (31)-(33) take the compact form of the Liouville-von Neumann equation 

ihd t p{k, t) = [Ho(fc) + S(fc, t) - tl(k, t), p(k, t)} (37) 

where [..,..] indicates the antisymmetrized matrix product. 

V. DISCUSSION 

This Section is devoted to the demonstration that Eqs. (31)-(33) provide the theoretical frame for describing different 
scenarios of current interest which are related to excitation with arbitrary light polarization and spin relaxation. These 
scenarios will be (i) the coupled two-level systems considered by San Miguel et al. in order to treat the polarization 
dynamics in VCELs, (ii) the quasi-equilibrium limit leading to a generalized inhomogeneous equation for the interband 
polarization matrix and (for the low-density limit) to the multiband Wannier-exciton equation, (iii) spin relaxation 
of optically oriented electrons which at the same time are due to the DP and BAP mechanisms, and (iv) the circular 
photogalvanic effect with interband excitation. 

(i) Following San Miguel et al. 10 in their phenomenological treatment of polarization dynamics in VCSELs one has 
to consider only the two-level systems consisting of heavy hole (m v = ±3/2) and electron states (m c = ±1/2) which 
are dipole coupled by circularly polarized light. Thus Eqs. (31)-(33) are simplified by neglecting the light-hole states 
and the band mixing. For circular polarization the matrix of Rabi frequencies (without renormalization) takes the form 
^m c m„ ~ ^m c ±i/2^m„±3/2- The diagonal terms of Eqs. (31) and (32) can be combined to write the equation-of-motion 

for the population inversion D = h{p\i + — (2 — p\s — P_3_s)} as 

22 22 22 22 

d t D = -— Z? + |lm(Pi3ffi 3 +P_i_30* 1 3) 

Tj 1 h 2222 22 2 - 2 
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n* mcmv (k,t) = E(t)-d% cmv (k)+ ]T v^ m ,_ mvm ,(k- 



(38) 



where radiative decay is considered phcnomenologically by the first term on the right hand side. This equation is 
(except for cavity effects, not considered here, and after adapting the notation) identical with Eq. (2.5) of 10 . Likewise 
the equation-of-motion for the spin polarization d = o{p|i + /°t a — P < Li_i ~ 3_a} is obtained in the form 

2 2 22 22 22 

d t d= -— rf+|lm(Pi3ffi 3 -P_i_sfi* i s ) 

T|| h 22 2 2 2 2 2 ~2' 

+ 4fs C 1 _l j O e _H -p e i_lE= n -^_3^3_3 +S"_3 3 /_3 3 ) • (39) 
III \ 2 2 22 22 22 2222 22 22/ 

In contrast to Eq. (38), additional terms containing off-diagonal elements of the density and self-energy matrices occur 
on the right hand side. When identifying these terms as ——d, describing spin relaxation 33 , one recovers Eq. (2.6) 
of 10 . Finally Eq. (33) can be formulated as 

d t P H±i = -^ +iuJ )p H±i + ^ H±i (D±d) 

-j(^^ ± |-^|^ i±i ) (40) 
where %uj — e c , i {k) — e v , 3 (k) and a phcnomcnological damping of the interband polarization was added. The last term 

± 2 ± 2 

(having a similar structure as that of Eq. (39)) can be understood as a contribution to this damping and combined 
with the first term to obtain Eq. (2.4) of 10 . Thus the phenomenological Maxwell-Bloch equations (especially those 
of Ref. 10 except for cavity effects) result as a special case of Eqs. (31)-(33) and give them a microscopic justification 
including an interpretation of the spin- relaxation time t s , which will become more transparent in the following. 

(ii) Equations (31)-(33) in combination with Eqs. (27)-(29) involve different time-scales in the evolution of the density 
matrices for electrons and holes. Scattering due to terms neglected in the HF truncation lead on the sub-ps time-scale 
to a loss of coherence concomitant with separate thermal equilibrium of electrons and holes (of different spin) that 
decays by spin relaxation and electron-hole recombination. This quasi-equilibrium can be described by replacing 

Pm c m c (k, t) - f° mc (k)5 m ^ c , ft v rn v *) ^ (1 " .C„ (k))S mv rn v (41) 

where /„(&), (&) are the Fermi-Dirac distribution functions for conduction and valence band states, with quasi- 
Fermi levels depending on m c ,m v . For this case the complex conjugate of Eq. (33) may be reformulated by making 
partial use of the matrix notation (introduced at the end of section IV): The matrix elements P mc m„ form the matrix 
of the interband polarization P(fe, t), (e^ c (fc) — s^ nv (k))P mc m v (k,t) can be replaced by [Ho(fc), P(k, i)] TOcTOt) , and the 
terms containing self energies by [S(fc, t), P(fc, i)]m c m„ • Thus, Eq. (33) takes the form 

ihd t P mcmv (k,t) - [Ho(fc) + £(M),P(M)] mcm „ = {f c mc {k) - f v mv {kWm c m v {k,t). (42) 

As in Eq. (40) (to which it is related) one could add a phenomenological damping term. Equation (42) is the general- 
ization of the RPA equation for the interband polarization (see Ref. 13 ) to the multiband case. 

It is interesting here to simplify to the case without band-mixing, for which the electron-electron and hole-hole 
interactions become diagonal in the spin indices in contrast with the electron-hole exchange interaction. The latter 
contributes to the off-diagonal elements of the electron self-energy matrix S c (/c) for nonvanishing hole concentration, 
thus indicating its relation to the BAP mechanism of spin relaxation. 22 Further specialization to the low-density limit 
leads to the multiband generalization of the Wannier-exciton Schroedinger equation. 

(iii) Spin relaxation of electrons due to spin-orbit coupling combined with momentum scattering is described usually 
starting from an equation-of-motion for the electron spin-density matrix p e (k, t). 5 ' 7 Such an equation is obtained from 
Eq. (31) by applying matrix notation as in Section IV but now only for the 2x2 electron block and identifying the 
first term on the r.h.s. as the generation matrix G(k, t) (see also Ref. 34 ) 

dtP e (k, t) - ^[HS(fc) + S c (fc), p e (k, t)} + W(k, &)(p e {k, t) - p e (fc', i)) - G(fc, t) , (43) 

fe' 

where we have added the last term on the l.h.s. to account for momentum scattering with phonons and (nonmagnetic) 
impurities. Equation (43) is a generalized form of equations used in the context of spin relaxation due to the DP 
mechanism in that it combines the single-particle spin- relaxation mechanism due to spin-orbit coupling (DP) with the 
many-body mechanisms considered in the electron self energy S c , in particular the one caused by the electron-hole 
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exchange interaction (BAP). To the best of our knowledge such a unified description of both mechanisms does not 
yet exist in the literature. A similar equation following from Eq. (32) for the 4x4 hole-density matrix p h (k, t) could 
be used to calculate the hole-spin relaxation. 

(iv) Reported measurements of the photogalvanic effects (PGE) in QW structures have been performed with mid- 
and far-infrared light on p- and n-doped samples. 16 ' 18 The nonlinear excitation with elliptic polarization creates non- 
equilibrium populations of different spin states by intersubband transitions which in combination with scattering 
processes result in a stationary current density described by 

j a = X af 3jE fj (w)£7 7 (-w) . (44) 

It is ruled by a third rank tensor (similar to \^ for second harmonic generation) that has nonvanishing components 
only in (so-called gyrotropic) systems lacking inversion symmetry. It can be decomposed into contributions from 
linear and circular polarization. An alternative formulation of this stationary current density is given by 34 

j a = eTr(pv a ) (45) 

where v a is a component of the velocity operator. As discussed in the literature the evaluation of Eq. (45) requires 
to consider the off-diagonal elements of p and v a , where p contains the nonlinear dependence on the electric field 
amplitude of the exciting light. The theory developed here does not directly apply to this situation but can be adapted 
to account for nonlinear intersubband excitation. On the other hand investigations of the CPGE with optical (or 
valence to conduction band) excitation, to which the concept presented here applies, are conceivable. An aspect of 
particular interest of PGE measurements is the dependence of the saturation behavior on the light polarization, which 
as an alternative to time-resolved experiments provides the possibility to detect the spin-relaxation time. Future work 
has to show how this dichroism of the PGE saturation can be described in the frame of the concept presented here. 
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